The Effect of Hot Gas in WMAP's First Year Data 
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i!B ■ ABSTRACT 

I ' By cross-correlating templates constructed from the 2 Micron All Sky Survey (2MASS) Ex- 

O I tended Source (XSC) catalogue with WMAP's first year data, we search for the thermal Sunyaev- 

-4— » ' Zel'dovich signature induced by hot gas in the local Universe. Assuming that galaxies trace the 

jrt ! distribution of hot gas, we select regions on the sky with the largest projected density of galaxies. 

Under conservative assumptions on the amplitude of foreground residuals, we find a temperature 
. ^_, , decrement of -35 ± 7 fiK (~ Str detection level, the highest reported so far) in the ~ 26 square 

^ ' degrees of the sky containing the largest number of galaxies per solid angle. We show that most 

JH I of the reported signal is caused by known galaxy clusters which, when convolved with the average 

beam of the WMAP W band channel, subtend a typical angular size of 20-30 arcmins. Finally, 
after removing from our analyses all pixels associated with known optical and X-ray galaxy clus- 
ters, we still find a tSZ decrement of -96 ± 37 fj.K in pixels subtending about ~ 0.8 square degrees 
on the sky. Most of this signal is coming from five different cluster candidates in the Zone of 
Avoidance (ZoA), present in the Clusters In the ZoA (CIZA) catalogue. We found no evidence 
that structures less bound than clusters contribute to the tSZ signal present in the WMAP data. 

Subject headings: cosmic microwave background — large-scale structure of universe — galaxies: clusters: 
general 

1. Introduction the CMB temperature field with unprecedented 

accuracy (Spergel et al. 2003). This tempera- 

The study of the Cosmic Microwave Back- ^^^.g ggj^ is composed of primordial anisotropics, 

ground (CMB) has become a powerful cosmologi- generated at the Last Scattering Surface, and sec- 

cal tool with applications in various astrophysical ondary fluctuations, which arise as the CMB pho- 

scenarios. Recently, the WMAP team has de- ^^^g ^^^^^^1 ^ ^^^ observer. Among these sec- 
termined the main cosmological parameters from 



ondary anisotropies, we shall concentrate on the 
so called thermal Sunyaev-Zel'dovich effect (here- 
after tSZ, Sunyaev & Zeldovich (1980)), associated 
with the distortion of the black body spectrum of 
the CMB photons due to Compton scattering on 
fast moving thermal electrons. This spectral dis- 
tortion is independent of redshift and proportional 
to the integrated electron pressure along the line 
of sight. It has been detected in the direction 
of several galaxy clusters (Carlstrom et al. 2002). 
This makes the tSZ effect an useful tool for the 
detection of ionized hot gas. 

Recently, Fukugita & Peebles (2004) have ar- 
gued that about 90% of all baryons are in the 
form of intergalactic plasma, not dense and hot 
enough to be detectable as X-ray sources except 
in clusters and groups of galaxies. At low redsh- 
fits, 30% of all baryons could be in the form of 
Ly-a absorbers (Pcnton, Stockc & ShuU 2004), 
but a large fraction of them has not yet been ac- 
counted for obscrvationally. Most models of struc- 
ture formation predict baryons to be located in fil- 
aments and sheets, associated to galaxy overden- 
sities (Springel, White, & Hernsquist 2001). KuU 
& Bohringer (1999) show evidence of filamentary 
X-ray emission at the core of the Shapley Super- 
cluster, whereas Zappacosta et al. (2002) report a 
detection of diffuse X-ray emission by warm gas 
(T ~ lO^K) associated with an overdense galaxy 
region. The aim of this letter is to use the tSZ 
effect to detect directly the diffuse warm baryon 
component in the local Universe. The analyses 
of the first year WMAP data have indicated that 
a small contribution due to the tSZ induced by 
clusters is present in the data (Bennett et al. 
(2003), Hernandez-Montcagudo & Rubiiio-Martin 
(2004)). Fosalba, Gaztaiiaga & Castander (2003) 
and Afshordi, Loh, & Strauss (2004) claimed a 
tSZ detection at small angular scales but did not 
clarify the nature of the astrophysical sources as- 
sociated with it. Myers et al. (2004) analysed 
the angular extension of the tSZ signal in the 
W band of WMAP by examining the CMB map 
in the direction of the largest galaxy groups and 
clusters found in the catalogues of Abell, Corwin 
& Olowin (1989) (hereafter AGO), Dalton et al. 
(1997) (hereafter APM) and the 2 Micron AU Sky 
Survey Extended Source Gatalogue, (Jarrett et al. 
2003). Their analyses showed evidence for diffuse 
tSZ emission up to an angular scale of '^ 1°, which 



implied a baryon fraction larger than the WMAP 
estimate. 

If baryons are distributed like the dark mat- 
ter on scales comparable to the virial radius of 
galaxies (Fukugita & Peebles 2004), the projected 
density of galaxies on the sky will correlate with 
the tSZ signal, independently of whether galaxies 
form clusters, groups or filaments. In this letter, 
we shall carry out a pixcl-to-pixel comparison of 
WMAP W band data with templates constructed 
from the 2MASS galaxy catalogue, to test the dis- 
tribution of hot gas in the local Universe. In Sec. 2 
we describe the pixel to pixel comparison method 
employed to estimate the tSZ signal and the data 
used in the analysis. In Sec. 3 we present and dis- 
cuss our results, and we conclude in Sec. 4. 

2. Method and Data Sets. 

The brightness temperature measured by a 
GMB experiment is the sum of different com- 
ponents: cosmological Tcmb, tSZ, instrumental 
noise N and foreground residuals F. If the tSZ 
signal is well traced by a known spatial tem- 
plate, denoted here as M and built for exam- 
ple from a galaxy catalogue, the total anisotropy 
at a fixed position on the sky can be modelled as 
T = Tcmb + a-M/{M)+N + F, where a measures 
the amplitude of the template induced signal and 
(M) denotes the spatial average of the template. 
If all other components have zero mean and well 
known correlation functions, then it is possible to 
use a pixel to pixel comparison to estimate a (see 
Hernandez-Monteagudo & Rubifio-Martin (2004) 
for details). If C denotes the correlation matrix 
of the GMB and noise components (foregrounds 
residuals will be discussed later) then the estimate 
of a and its statistical error are 
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Since our galaxy template M will be positive by 
construction, this equation demands that GMB 
and noise fields have zero mean. We impose the 
average of all pixels outside the KpO mask to be 
zero. We checked that our results were insensi- 
tive to this requirement by carrying out a similar 
analysis usig the more conservative Kp2 mask. In 
all cases our results changed by less than a few 
percent. Notice that our method requires the in- 



version of the correlation matrix, a computation- 
ally expensive procedure. To speed up the process 
we will carry out the analysis in pixel subsets as 
described below. We checked with Monte Carlo 
simulations that (Tq, is an unbiased estimator of 
the error. 

We centered our analyses on WMAP W band, 
which has the highest angular resolution. In this 
band, the instrumental noise shows almost no spa- 
tial correlation and its position dependent am- 
plitude is determined by the number of observa- 
tions^ (Jarosik et al. 2003). The tSZ template was 
built from the 2 Micron AU Sky Survey (2MASS) 
Extended Source Catalogue (Jarrett et al. 2003). 
This catalogue contains about 1.5 million galaxies, 
extending up to redshift z ^ 0.1 (400 Mpc), de- 
tected in the near- infrared (J, H and Kg bands). 
In order to build a tSZ template (M), aU 2MASS 
galaxies were projected onto the sphere using the 
HEALPix^ pixehzation (Gorski, Hivon & Wan- 
delt 1999) with the same resolution as the CMB 
data. The amplitude in every pixel was made pro- 
portional to the number of galaxies. The tem- 
plate was then convolved with the window func- 
tion of the noise weighted average beam of the 
four Difference Assemblies (DA's) corresponding 
to the WMAP W band (clean map), and multi- 
plied by the KpO mask, like the CMB data. Since 
our method requires the CMB and noise compo- 
nents to have zero mean, we substracted the map 
average outside the KpO mask. 

Pixels were sorted in sets of size Npi^, de- 
noted as M^, where the superscript index [i in- 
dicates galaxy density, in such a way that low 
(3 corresponds to higher projected galaxy density. 
For patches of Npix =2048, the average projected 
galaxy density for (3 =1 was ^ 420 galaxies per 
square degree, whereas for (3 =500 the density 
dropped to ~ 44 galaxies per square degree, which 
roughly coincides with the average projected den- 
sity outside the KpO mask. Npix ranged from 
64 to 2048. The pixel to pixel comparison was 
then performed on each of these subsets: we com- 
pared all pixels in M^ to their counterparts in the 
CMB map. Our working hypothesis is that galax- 
ies are fair tracers of the gas density and within 



each set electron temperature is similar, that is, in 
those pixels galaxies trace the electron pressure. 
Within each subset M'^, the galaxy density re- 
mained roughly constant, so our method returned 
a as a weighted mean of the measured temperature 
in those pixels. If no tSZ signal is present, then a 
will scatter around zero, as a consequence of CMB 
and noise being random fields of zero mean. 

3. Results and Discussion. 

Fig. (1) summarizes our main results: Fig. (la) 
shows the estimated a's for the sets having the 
highest projected density of galaxies. In abscis- 
sas we give the set index (/3). Crosses, filled cir- 
cles, triangles and diamonds correspond to Npix — 
256, 512, 1024 and 2048, respectively. Symbols 
are slightly shifted for clarity, error bars denote 
Ict confidence levels. At WMAP instrumental fre- 
quencies, the tSZ effect causes temperature decre- 
ments, i.e., template and CMB data should an- 
ticorrelate giving negative a's, as found. Sets 1 
and 2 of size Nr,i 
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pix 



iThe WMAP data was 

http: //lambda, gsfc.nasa. gov. 
^ http: //www. eso. org/science/healpix/ 



downloaded 



from 



: 256 correspond to the first 
512 and similarly for all other sizes 
and sets. Consistently, a of a larger set is always 
bracketed by the a's measured from its subsets. 
The largest signal comes from the densest 256 pix- 
els but the highest statistical level of significance is 
achieved for Npi^ — 2048 (a — —35 /iK at the 4.9(T 
detection level), since the error bars shrink due to 
the higher number of pixels contributing with tSZ 
signal. In Fig. (lb) the same data sets are plotted 
versus the average projected galaxy density, i.e., 
the average number of galaxies per square degree 
within each subset, as seen by the W band beam. 
The symbol coding is identical: first sets of Npix — 
256, 512 show a projected galaxy density as high 
as 420-380 galaxies per square degree, whereas the 
6th set for Npix ~ 2048 contains around 140 galax- 
ies per square degree. 

In order to evaluate the significance of the pre- 
vious results, we repeat the analysis for pixels of 
intermediate and low projected density of galaxies. 
In Fig. (Ic) diamonds correspond to sets of size 
Npix = 2048 with (3 E [1, 30], whereas filled circles 
correspond to indices /3 G [501,530]. The shaded 
area limits the la error bar for diamonds, that 
is about a factor of 1.2 bigger than for filled cir- 
cles. While the latter scatter around zero with the 
expected dispersion, diamonds are clearly biased 



towards negative values: besides the first patch 
(a = —35 ±7 fJ-K), there are seven other sets above 
the 2(7 level, one of them at 3cr. 

We also tested the consistency of our results 
with respect to frequency. In Fig. (Id) we esti- 
mate a by cross-correlating the densest pixel set 
{Npi^ = 2048, /? == 1) with all WMAP bands: 
K (23GHz), Ka (33GHz), Q (41GHz), V(61GHz) 
and W(93GHz). Diamonds give the obtained a's 
from raw maps, whereas triangles refer to analy- 
ses performed on foreground cleaned maps avail- 
able in the LAMBDA site. Note that the WMAP 
team only provided clean maps for the three high- 
est frequency channels. In all maps, the signal is 
compatible with being due to tSZ. A more quanti- 
tative comparison is not straighforward since the 
maps have different angular resolution and galac- 
tic contamination. When compared in pairs, all 
a's are within Icr of each other and they are all 
within 2cr of the expected frequency dependence 
of the tSZ effect, whose best fit for the cleaned Q, 
V and W band maps is plotted as a solid line. 

We measure the extent of the tSZ regions by 
rotating the template around the z-axis (perpen- 
dicular to the galactic plane) and cross-correlating 
the densest 2048 pixels. In Fig. (2) squares give 
a versus the average angular displacement of the 
pixel set. Error bars are again la. The solid line 
represents the gaussian approximation of the W- 
band beam. The size of the tSZ sources is typically 
20-30 arcmins, slightly bigger than the beam but 
remarkably smaller than the values obtained by 
Myers et al. (2004) from rich AGO clusters. Actu- 
ally, when one studies the angular distribution of 
the 2048 densest pixels, one finds that they are as- 
sociated in groups of typically 3-4 members, and 
that the groups are uniformly distributed on the 
sky. 

By comparing, in Fig. (Id), results from the raw 
and cleaned maps built from the W band, one can 
conclude that foregrounds will have little impact 
on our results. We made a more detailed study 
using Monte Carlo simulations: we performed 100 
simulations of the GMB and W-band noise compo- 
nents and added and removed a foreground resid- 
ual template. As a conservative model for fore- 
ground residuals left after cleaning the W-band, 
we took the sum of the dust, free-free and syn- 
chrotron emission maps for the W band released 
by the WMAP team. In each simulation we (i) 



added and (ii) substracted the whole foreground 
template. In Fig. (2) we plot the average a's for 
these simulations, after adding (dotted line) and 
substracting (dashed line) the residuals for several 
angular displacements. The shaded bands display 
the la dispersion areas. The errors practically 
equal the statistical estimates of eq. (1). To sum- 
marize, foregrounds do not significantly affect our 
estimated amplitude of the tSZ contribution. 

Finally, we removed from our galaxy template 
all those pixels that were associated with known 
galaxy clusters. We used the AGO and APM cat- 
alogues of optically selected clusters and the XBC 
(Ebeling et al. 2000), de Grandi (de Grandi et al. 
1999), NORAS (Bohringer et al. 2000), ROSAT- 
PSPG (Vikhlinin et al. 1998) and Vogues (Vogues 
et al. 1999) X-ray cluster catalogues. We excised 
from the analyses all pixels lying within a virial 
radius of the cluster center (taken to be ten times 
the core radius). For clusters without measured 
core radius but with known redshift, we assumed 
a virial radius of 1.7 Mpc, and removed all pixels 
within that distance. For the rest (the majority), 
we conservatively removed all pixels within a circle 
of 30 arcmin from the cluster centre. Out of the 
2048 pixels ioi (3 ~ 1, 1681 were eliminated. For 
the patches with (3 — 2 — 30, i.e., the next ^60,000 
densest pixels in Fig.(lc), a large fraction of them 
were also associated to known clusters and were 
eliminated with the excising. In Fig. (3) we show 
the cross correlation of these remaining pixels out- 
side known clusters with the clean W band map. 
Note that here patch sortening has been regener- 
ated using surviving pixels. For the densest sets of 
Npix = 64 (filled circles), 128 (triangles) and 256 
(diamonds), we still find evidence of tSZ, but at 
much lower level of significance. For the densest 64 
pixels, subtending ~ 0.8 square degrees on the sky, 
we obtain a = — 96±37 /iK, at ~ 2.6cr significance 
level. The signal gets dilluted rapidly as more pix- 
els are included in the analysis: a ~ —50 ± 27 /iK 
and a = -30 ± 18 fiK for Npi^ = 128, 256, re- 
spectively. Out of the 64 densest pixels, 54 pixels 
are in the ZoA, and 45 of them coincide with five 
different cluster candidates in the GIZA (Ebeling, 
Mulfis, & Tully 2002) catalogue, (Ebeling, private 
communication). The remaining group of pixels 
are not associated to any known galaxy cluster. 
In Fig. (4) we plot the location of those pixels in 
the sky. The shaded area correspond to the KpO 



mask, that also remove many point sources outside 
the galactic plane (dark grey dots). The 64 pixels 
are ploted as big white circles for convenience. 

4. Conclusions 

Under the assumption that galaxies trace hot 
gas, we have used the 2MASS galaxy catalogue 
to search for tSZ signal present in WMAP data. 
In ^ 26 square degrees on the sky we have found 
a contribution of average amplitude -35 ± 7 /iK, 
spectrally compatible with tSZ, and mostly gen- 
erated by AGO clusters of galaxies. Our study, 
based on a pixel-to-pixel comparison, reaches the 
highest sensitivity level reported so far. Compared 
with methods based on power spectrum or corre- 
lation function analysis, our method gives a larger 
level of significance since we restrict the analysis to 
the regions of the sky where the tSZ contribution 
is expected to be the largest. 

We have found that the typical angular exten- 
sion of this signal is somewhere between 20-30 ar- 
cmins. Furthermore, once all known clusters of 
galaxies are excised from the analysis, we are left 
with ^ 0.8 square degrees with an average ampli- 
tude of q; = —96 ± 37 /iK: those pixels fall mostly 
in the ZoA and after performing our analyses we 
found that 45 of them are associated to five dif- 
ferent galaxy clusters in the CIZA catalogue. We 
have found no conclusive evidence that, in the vol- 
ume probed by 2MASS, structures less bound than 
clusters contribute to the tSZ signal present in the 
WMAP data. 
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Fig. 1. — (a) tSZ detection in pixel sets with the highest projected galaxy density. Crosses, filled circles, 
triangles and diamonds correspond to sets with Npix — 256, 512, 1024, 2048, respectively. Bars represent la- 
errors, (b) tSZ signal versus projected galaxy density. The symbol coding is the same as in panel (a), (c) 
tSZ signal for the very dense (M^ to M^°, diamonds) and less dense (M^°^ to M^'^", filled circles) sets. 
Bars and band are Icr errors, (d) tSZ for the different WMAP bands estimated on the densest set of 2048 
pixels. The solid line represents the expected frequency dependence fitted to the clean Q, V and W band 
maps (triangles). Estimates of a obtained from raw (i.e., non cleaned) maps are displayed by diamonds. For 
bands with clean and raw maps, triangles and diamonds are slightly shifted for better display. 




Fig. 2. — Rotation test showing the angular ex- 
tent of the reported tSZ sources: a is computed 
for different angular displacements of the galaxy 
template with respect to the cleaned W band map, 
for the densest 2048 pixels (squares, error bars at 
la). This analysis is repeated on 100 MC simula- 
tions including CMB, W-band noise plus (dotted 
line) and minus (dashed line) our foreground tem- 
plate. Shaded areas denote again Icr confidence 
level. The solid lines displays the gaussian ap- 
proximation to the W band beam. 
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Fig. 3. — Persistence of tSZ signal once all 
known clusters have been excised from the anal- 
ysis. Filled circles, triangles and diamonds cor- 
respond to patches of Npi^ — 64, 128 and 256, 
respectively. 




Fig. 4. — White spots show the position of the 64 
densest pixels not excised after removing known 
galaxy clusters. These pixels contribute with a 
tSZ signal of amplitude of a = —96 ± 37 fiK. The 
dark grey area is masked out by the KpO mask, 
which covers the Galaxy and the brightest radio 
sources. The graticule is 20° x 20°. 
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